Kondo resonance for orbitally degenerate systems 
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Formation of the Kondo state in general two-band Anderson model has been investigated within 
the numerical renormalization group (NRG) calculations. The Abrikosov-Suhl resonance is essen- 
tially asymmetric for the model with one electron per impurity (quarter filling case) in contrast with 
the one-band case. An external magnetic (pseudo-magnetic) field breaking spin (orbital) degener- 
acy leads to asymmetric splitting and essential broadening of the many-body resonance. Unlike the 
standard Anderson model, the "spin up" Kondo peak is pinned against the Fermi level, but not 
suppressed by magnetic field. 
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The Kondo problem is one of most fascinating and im- 
portant issues in condensed matter theory. It was orig- 
inally formulated for explanation of the resistivity mini- 
mum in metallic alloys |l| due to the scattering of con- 
duction electrons by a magnetic impurity, and later gen- 
eralized to various cases. The Kondo effect turns out 
to be a key phenomenon of the heavy fermion behavior 
, anomalous electronic properties of metallic glasses at 
low temperatures Q , quantum dots , and many other 
correlated electron problems. One of the main results of 
this theory is the formation of a resonance with small 
energy scale (so-called, the Kondo temperature Tk) near 
the Fermi energy due to the scattering of conduction elec- 
trons by local quantum systems with internal degrees of 
freedom. Originally this resonance (usually called the 
Kondo, or Abrikosov-Suhl resonance; for a review, see 
Ref. |2) could be experimentally investigated only indi- 
rectly, through the temperature dependence of thermo- 
dynamic and transport properties of metals. However, 
novel tunneling spectroscopy, in particular scanning tun- 
neling microscopy (STM) , is now able to directly visualize 
the Kondo resonance 0,0,0- At, the same time, the- 
oretical investigation of the spectral density for Kondo 
systems is a much more complicated problem than the 
calculation of thermodynamic properties. Indeed, in the 
latter case exact analytical results can be derived by the 
Bethe Ansatz H], whereas the electron density of states 
is considered mainly by some approx imate methods or 
numerically, see, e.g., Refs. II El 111 El El El El El 
The splitting of the Kondo resonance by external mag- 
netic field was investigated by Bethe Ansatz for the s-d 
exchange model E3- These results can be used to verify 
different approximate schemes demonstrating some dif- 
ficulties, for example, with the well-known non-crossing 
approximation (NCA) (lfij . 

Generally speaking, we have now a complete and satis- 
factory theory of the Kondo resonance for the prototype 
case of purely spin scattering. On the other hand, the 
information about the systems with orbital degrees of 



freedom is still insufficient. A so-called "orbital Kondo 
resonance" has been considered theoretically, for atomic 
two-level systems in metallic glasses, for quadrupolar de- 
grees of freedom in some uranium-based compounds H , 
for high-temperature superconductors [l8l |. and for dou- 
ble quantum dot systems |l9j . Recently the phase di- 
agram of the Anderson model with orbital degrees of 
freedom has been investigated b y th e numerical renor- 
malization group (NRG) method |2fj. Spin (pseudospin) 
susceptibility for the double quantum dot model has been 
investigated by this technique in Ref. El However, elec- 
tron spectral density has not been calculated. It is worth- 
while to note that investigation of dynamical properties 
for effective impurity models is of the crucial importance 
for the dynamical mean-field theory (DMFT) approach, 
in particular, to describe the metal-insulator transition 
and related phenomena |2l| . For one-band case the NRG 
method was applied to DMFT problem in Ref. |H 

Direct observation of the orbital Kondo resonance on 
Cr(OOl) surface by the STM measurements H , as well as 
a relevance for multiple quantum dot systems |lj| , makes 
the issue about the shape of the Kondo resonance and 
about effects of "pseudomagnetic" field which breaks the 
orbital degeneracy especially actual. Here we investigate 
the problem of orbital Kondo resonance by the NRG ap- 
proach [il |. We shall demonstrate that in the "orbital" 
Kondo case the resonance has an essential asymmetry 
with respect to the Fermi energy, in qualitative agree- 
ment with the experimental observations 0. 

We start from the two-band Anderson impurity model 
with the spin and orbital rotationally invariant Hamilto- 
nian El 
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where 



H imp - ^2 (ff - "2") /«L/< 



U + J 



Here a, 6 = 1,2 and cr =|, j arc orbital and spin indices, 
correspondingly, (cka<r) denote the creation (anni- 
hilation) operators for a-orbital states with spin a and 
energy £k (we take the rectangle band with half width 
D = 2), faaifao) those for impurity states of a-orbital 
with spin a and energy £f, n acr = /J^/aa, h is magnetic 
field. Since the spin and orbital degrees of freedom are 
symmetric in the Hamiltonian, h may be a pseudomag- 
netic field (e.g., for the orbital Kondo effect on Cr(OOl) 
surface where the potential of the atomic step edge breaks 
the exact degeneracy between d xz and d yz states Q). 
The Coulomb interaction and exchange parameter at the 
impurity site are U and J, and two electron subsystems 
are coupled via a hybridization parameter V. Note that 
in solids the orbital moment conservation is not an ex- 
act property and therefore some additional terms in the 
Hamiltonian may appear |25l |. However, we omit them 
since the problem turns out to be numerically very cum- 
bersome even for the rotationally invariant Hamiltonian. 

To calculate the spectral properties of impurity we use 
the NRG technique which is described in details in Refs. 
0, El EE EE Here we emphasize the new aspects for 
multi-orbital Anderson model. 

As usually, we start from the solution of the isolated 
impurity problem. As an initial step of the RG pro- 
cedure, we add the first conduction electron site, diag- 
onalize the Hamiltonian matrix for this Hubert space, 
and thus obtain new eigenstates. Then such a procedure 
should be repeated until reaching a fixed point. Dimen- 
sion of the Hilbert space within one NRG-iteration in- 
creases by factor 16 instead of 4 for the one-band case. By 
using an appropriate symmetry of the problem we were 
able to find the whole spectrum for the Hilbert space 
with the dimensionality of 48000 states, which give s a 
possibility to keep at each step about 3000 states |27j . 

Due to the NRG discretization scheme, the electron 
spectral function 



p f (uj) = ImG(w + iO) 
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is given by a set of ^-functional peaks at the frequen- 
cies oj n . Standard NRG practice consists of the Gaussian 
broadening of the spectral function on a logarithmic scale 
|llj . Since the point u> = plays special role in such a 
scheme, we used more conventional Gaussian broadening 
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smearing being changed depending on the iteration num- 
ber L, namely, &l+2 = &l/A (where the NRG cutoff pa- 
rameter A = 2 have been used). 
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FIG. 1: The density of states for the half-filled case, et = 
—3, V = 0.2, U = 2, J — 0, h = (solid) and symmetric split- 
ting (only the spin-up DOS) in the magnetic field h — 0.006 
(dot), h — 0.01 (dash) and h = 0.02 (dash-dot). Occupation 
number per orbital and spin are shown in the left insert (total 
nf — 2). Density of states at the Fermi level according to the 
Friedel sum rule © is pf (0) = 10.13. Right insert shows the 
DOS in the vicinity of the Fermi level. 

The Kondo regime corresponds to the case where 
| £ f| ^ 2r [12{. r = irV 2 p being the one-particle reso- 
nance width for the localized electrons, p is the bare DOS 
of conduction electrons at Ep. The density of states for 
the half-filled case (n/ = 2) and for almost quarter-filled 
case [n t ~ 1) are shown in Figs. 1 and 2, respectively. 
Although the Hund interaction J was estimated to have 
considerable value 0, this quantity is irrelevant for the 
quarter-filled case; our RG calculations confirmed that 
J does not influence essentially the results (see insert in 
Fig. 2) . The main difference of the NRG results obtained 
here with those for the nondegenerate Anderson model is 
that the Kondo peak is not centered at the Fermi energy 
Ep = 0. An explanation of this deviation is given by the 
Friedel sum rule for the phase shifts r]f. 2(2l+l)r]i/n = 71/ 
E^ ,21 + 1 being the number of orbital channels and n / 
the total number of localized electrons. It is important 
that due to locality of the self-energy the value of p/(0) 
does not change in comparison with the non-interacting 
model (U = 0) and is equal to 
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where N is the degeneracy factor |2jj . For the standard 
SU(2) Kondo model for S = 1/2, as well as for a degen- 
erate half-filled model, the phase shift at Ep is close to 
7r/2 which means the strongest possible scattering at the 
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top of the resonance peak. In this case the asymmetry 
of the Kondo resonance with respect to Ep should be 
very weak even for the nonsymmetric Anderson model, 
in agreement with recent computational results [Hi ] , and 
with our Fig. 1. Thus the large value of electronic effec- 
tive mass and linear specific heat is owing to renormal- 
ization of the residue of the electron Green's function, 
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On the other hand, for N > 2 and n/ « 1 the top of the 
f-electron peak shifts above Ep In particular, in the 
SU(N) Anderson model (or in the equivalent Coqublin- 
Schrieffer model) with infinitely large N 0, 0| the high 
peak lies completely above the Fermi level, being shifted 
by the value of the order of Tk (note that in the limit 
N — > oo this is just a delta-like peak). NRG calcula- 
tions for SU{N) Anderson model [llj give qualitatively 
similar results. However, the values of p/(0) calculated 
in Ref. 01 turn out to be by about 30% smaller than 
given by Eq.|j3J), probably, due to insufficient accuracy. 
Our results for double-degenerate model give p/(0) in the 
perfect agreement with the Friedel sum rule (see captions 
to Fig. 1). 




FIG. 2: The density of states for the almost quarter-filled 
case for et = —1, U = 2, J = 0, V = 0.35. Occupation num- 
ber per orbital and spin is 0.258 (total nj = 1.032). Insert 
shows the DOS in the vicinity of the Fermi level; the dash 
line corresponds to the case of non-zero exchange parameter 
J = 0.2. 

The NRG calculations with the external magnetic field 
h were also performed. Since we are interested only in 
a small energy region, shortcomings of the NRG scheme, 
which occur in the presence of magnetic field at large ex- 
citation energies to |3dl| . are not important here. There 



are several qualitatively different regimes of magnetic 
splitting in the double-degenerate Anderson model: 

(i) a half-filled case where n/ = 2 due to electron- 
hole symmetry even in the presence of external magnetic 
field. The usual symmetrical splitting takes place, the 
Kondo peaks becomes low and broad (Fig. 1), similar to 
the behavior in the non-degenerate symmetric Anderson 
model [3. 

(ii) a nearly quarter-filled case: Uf f=a 1. The split- 
ting of the Kondo peak is asymmetric. One can see from 
Fig. 3 that the upper "spin down" peak becomes low 
and broad, as well as in the symmetric case. At the same 
time, the lower "spin up" peak is not suppressed (as in 
the standard Kondo model) , but tends to the Fermi level 
and becomes more high (however, the situation changes 
with decreasing V: the height and area of spin-up peak 
decreases considerably in strong magnetic field). The 
density of states at Ep remains high up to very strong 
fields, so that the partial f-occupation numbers depend 
strongly on h. For very strong fields the peak corre- 
sponding to spin down states is completely suppressed 
(jifai —>■ 0), and only orbital Kondo resonance between 
spin up states survives since Uf a \ — > 1/2, in contrast with 
the case (i) where n/ f — > 1. 

(hi) the intermediate valence regime (rif is essentially 
non-integer) . Instead of the three-peak structure charac- 
teristic for the Kondo regime, we have one peak which is 
split in magnetic field (Fig. 4). 
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FIG. 3: The effect of the magnetic field on the density of 
states for the same model parameters as in Fig. 2. The mag- 
netic field values are h — (solid), 0.006 (dash), 0.01 (dot), 
0.016 (dash-dot). The left insert displays the dependence of 
occupation numbers on the magnetic field; the right one shows 
the result of the approximation JSJ for gf = — 1, U — <x,J = 
0,V= 0.56, h = (solid) and h = 0.04 (dot), the smearing of 
the logarithm with S — 0.004 being introduced. The value of 
V is increased in comparison with the finite-?/ case to obtain 
a comparable value of the Kondo temperature. 

The NRG method can provide a detailed information 
about the Kondo resonance in the two-band Anderson 
model. For larger number of orbitals, numerical calcula- 
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FIG. 4: The effect of the magnetic field on the density of 
states for e f = -2, U = 2, J = 0, V = 0.2. The magnetic field 
values are h = (solid), 0.08 (dash), 0.32 (dot). The insert 
shows the spin-occupation number as function of magnetic 
fields. 



DMFT approach beyond the one-band case. Our version 
of the NRG can describe accurately the case where the 
Kondo peak is shifted from the Fermi energy, which is 
a generic case of multiband impurity model. The new 
features of the orbital-degenerate model are related with 
a fact, that the Kondo resonance is not suppressed by 
external magnetic (or pseudomagnctic) field, its splitting 
being essentially asymmetric. 
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tions become impractical for a reliable treatment of the 
Kondo problem. To understand qualitatively the behav- [5 
ior of the multiband Kondo system in magnetic field, 
we consider a simple degenerate Anderson model with [6 
U = oo where the peak lies also above the Fermi level 
[lfjj . The Hamiltonian reads 

r 1 [9 

U = zJ t ±4 m c km+52(£{-h m )f} n f m +V'%2 [cl m f m + h.cj 

km m km 

J [io 
[ii 

where f} n = |m)(0| are the Hubbard's AT-operators; the [12 
model (QJ corresponds to m — aa and h a -\ — h/2, h a i = P-^ 
—h/2. Usin g th e second-order perturbation theory for P^ 
X-operators |3jJ one can obtain (cf. also Refs. jd. Ilfl|) 

[16 

({fm\fm))E = (no + n m ) [E - Ef + h rn 

This very simple approximation gives reasonable quali- [20 
tative agreement with the numerically accurate NRG re- [21 
suits (see right insert in Fig. 3). It can be used for in- [22 
terpretations of the computational results. For example, 
according to Eq.JHJ), the spin up peak does not intersect 
the Fermi level with increasing magnetic field (at least, at 
not too small V). Indeed, the logarithmic divergence of 
the self-energy exists in our orbitally-degenerate model, [25 
the contribution from the transitions between the degen- 
erate states being not cut at h. [26 
To conclude, we have considered peculiarities of the ffl 
Kondo resonance in the orbital-degenerate case by a nu- 
merical renormalization group technique. A possibility to 
calculate the spectral properties for degenerate Anderson 
impurity model is demonstrated, which gives a chance to [30 
extend the applicability region of the NRG scheme in the 31 
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